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Abstract 

We examine whether the fact that QCD is chirally invariant at short dis- 
tances necessarily leads to the prediction that hadrons form approximate parity 
doublets, as has been recently put forward by Glozman and collaborators. We 
show that this is not the case, and we exhibit some of the pitfalls of trying to 
link the operator product expansion to the hadron spectrum. We illustrate our 
arguments within a model for scalar and pseudo-scalar mesons used recently by 
Shifman to argue for parity doubling. We find that, whatever the experimental 
situation may be, there is no theoretical basis for parity doubling based on the 
use of the operator product expansion. 
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1 Introduction 



In recent years there has been a renewed interest in the possibility that "chiral sym- 
metry restoration" takes place for the higher excited hadronic states of a given spin 
and isospin. For each spin and isospin, resonances occur with both parities, and if chi- 
ral symmetry, SU{rif)L x SU{nf)R, would be linearly realized these resonances would 
have to organize themselves in equal-mass pairs of opposite parity. Of course, it is well 
known that chiral symmetry is spontaneously broken, and that there is therefore no 
symmetry reason that such parity pairing should occur. This observation goes back to 
the work of Ref. £Q, and was recently re-emphasized in the language of effective field 
theory in Ref. |2j. 

However, one may ask whether there could be any other dynamical reason that 
approximately degenerate parity pairs occur in QCD. Loosely speaking, the intuition 
would seem to be that for highly excited resonances the physics of the vacuum (which 
breaks chiral symmetry) is less important. Such resonances would then couple less 
strongly to pions, and possibly organize themselves in approximate parity pairs. 

Here we examine one such dynamical scenario. 1 The basic observation is that, since 
chiral symmetry breaking is a non-perturbative phenomenon, and QCD is asymptot- 
ically free, chiral symmetry is restored at high energies, and, mutatis mutandis, for 
highly excited resonances [3]. To our knowledge, this has been most concretely ex- 
plored in Ref. [Sj in the context of the operator product expansion (OPE). 

In order to be precise, we begin with setting the stage for our discussion. First, in 
order to give an unambiguous meaning to the mass of a resonance, we will consider the 
limit in which N c , the number of colors, is taken to infinity. Otherwise highly excited 
resonances merge into the perturbative continuum, due to their finite widths, and their 
individual identity is lost. In contrast, in the limit of infinite N c , all resonances are 
stable, and spectral functions are (infinite) sums over delta functions. Second, we need 
a definition of what it means for chiral symmetry to be restored for highly excited 
resonances. We consider a tower of resonances of given spin, isospin and parity, with 
masses M$ n , and the corresponding tower with the same spin and isospin, but opposite 
parity, with masses Mp n . For definiteness, we will consider scalar and pseudo-scalar 
mesons, as was done in Ref. j3]. 2 The integer n labels the successive states in each 
tower. We will now assume that for each S state with label n, a closest P state can 
be identified, and for the sake of convenience we will use the same label n for this P 
state. 3 We now define as a measure of chiral symmetry the "spectral overlap" p(n) by 



The choice of denominator here deserves some comment. Obviously, we want p(n) to 
be dimensionless, but we do not want p(n) to go to zero for large n trivially. Therefore, 
the denominator is taken to be the smallest possible quantity with dimension of mass 

1 For a recent discussion of other possibilities, see Ref. 

2 Our arguments extend easily to the example of vector and axial-vector mesons. 

3 In other words, we assume that approximate parity doublets can be identified, at least for large 



6J 




M Pn - M s , 



(1) 



enough n. 
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made out of the masses with "resonance number" approximately equal to n. "Chiral 
symmetry restoration" for a particular spin and isospin takes place if p(n) goes to zero 
for large n. 

Another measure of "chiral symmetry restoration" considered in the literature |3 U\ 
is the "chiral asymmetry" 

, , M P - M s 
X{ } M Pn + M Sn 1 } 

This is not a good measure of chiral symmetry restoration, however, because of the 
choice of the denominator. We know that the resonance masses will increase with n, 
and thus this denominator will grow with n. That alone will make \ go to zero for 
large n, as long as the numerator grows less fast — in fact, typically one would envisage 
defining parity pairs such that the numerator stays finite for large n. Clearly, in order 
to avoid this problem, one should measure the difference in mass between potential 
parity partners relative to the overall mass density (per unit of n), as is done by the 
quantity p(ri). We will therefore not consider xi n ) i n this paper. 

The claim made in Ref. p. can now be expressed as follows: The slowest pattern of 
chiral symmetry restoration (for the towers of scalar and pseudo-scalar mesons) allowed 
by the asymptotic behavior of the OPE is given by 

\p(n)\ ~ - , n -»• oo . (3) 
n 

It was assumed that the asymptotic behavior for large n in both the S and P channels 
is Regge like, M Sn Pn ~ nA 2 . If Eq. (J3J) would be true, then clearly it would be 
appropriate to conclude that chiral symmetry is restored for highly excited resonances 
(in this channel). This result was inferred from a super convergence relation, 



(4) 



in which SM 2 = M| —Mp . Superconvergence relations of this type were also assumed 
to be valid in Refs. 0E|. ™ 

In this note, we show, within the context of the same assumptions (large N c and 
Regge behavior), that these claims do not follow from the OPE. 4 To clarify this point, 
we will, within the context of the same assumptions, exhibit a model of scalar and 
pseudo-scalar meson resonances which is consistent with what is known about the 
OPE, and for which p(n) approaches a non-zero constant for large n. 

This paper is organized as follows. We begin with reviewing the argument presented 
in Ref. j^j in the following section, and we point out why in general it is incorrect. We 
then remind the reader of the fact that large- N c QCD in two dimensions satisfies the 
assumptions, but violates the claims. It therefore is a good counter example, and our 
explicit calculation of p(n) for this model invalidates statements about this model in 

4 We have pointed out before that superconvergence relations such as Eq. (0J cannot be derived 
from the OPE |12| . However, in view of the recent claims which keep appearing in the literature, it 
seems worthwhile revisiting this issue once more. 
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Ref. In Sec. 4, we construct a model of scalar and pseudo-scalar resonances in four 
dimensions which further demonstrates that the claims of Ref. |5] do not have to hold. 
In our concluding section, we comment on the situation at finite N c . 



2 Scalar and pseudo-scalar mesons 

Following Ref. [5], we start with the scalar and pseudo-scalar two-point functions 

n 5 ,p(g 2 ) = i j dx e*** (0\T { J S ,p(z)4,p( )} 1°) > ( 5 ) 

with Js(x) = d(x)u(x) and Jp(x) = d(x)i^ 5 u(x), and consider a once-subtracted 
dispersion relation for one-half the difference of these two correlation functions: 

iW(g ) = — dt * — — — - — — — + An(o) g-v , (6) 

n Jo t(t - q z - le) -q 2 - te 

where the prime indicates that we do not include the pion pole inside the integral, 
instead writing it as a separate term. We work in the chiral limit, /o is the pion decay 
constant in that limit (approximately 87 MeV), and B = — (V'V'V/o- We note that 
there is no singularity in the integrals at t — because the spectral functions -Im Hs(t) 
and ^Im H'p{t) vanish at t = 0. The subtraction constant All, 

An(o) = l - (n s (o) - n' P (o)) = i6£? 2 L 8 (7) 

is fixed by chiral symmetry, where Lg is one of the order-p 4 couplings in the pion 
effective theory [TO] . 

Working at N c = oo, there is an infinite tower of infinitely narrow resonances in 
each channel, and we thus have that 

. oo 

-lmU s (t) = 2j2F 2 sAt-Ml), (8) 

n 

- x 

-Imn' P (t) = 2Y,F 2 P J{t-Ml n ) . 



Furthermore, we will assume Regge-like behavior asymptotically in n 0, i.e. 

M l,p n = nA 2 + • • • , n oo , (9) 
with A ~ 1 GeV, and we also require that 

*tA ~ A2M l,Pn . ^ ^ °° ' ( 10 ) 
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in order guarantee the correct parton- model behavior. 5 Substituting these expressions 
into Eq. (JBJ), one obtains 



u s Mq=-Q 2 ) = (ii) 

B 2 fi ^( Fl FJL 



Q 2 V \ M l (Q 2 + M l) M k (Q 2 + M k) 

which is finite, with the masses and residues of Eqs. © and |T0|) . 

Next, we define the mass splitting in a parity pair (for convenience we label both 
the S and P states of the parity pair with the same resonance number n, as already 
mentioned in the previous section) 

5M 2 n = Ml n -Ml n , (12) 

and we assume that SM 2 grows less rapidly than M 2 = (M| n + Mj> J /2 for n — > 
oo. Expanding Eq. (|TT|) in 5M 2 then yields, with Fg Pn = nA 2 Mg n P , where k is a 
proportionality constant fixed by the parton-model logarithm, 

This is to be compared with the leading term for large Q 2 in the OPE for this two-point 
function: _ 

IW(g 2 = -Q 2 ) ~ ^ , Q 2 - oo , (14) 

up to logarithmic corrections. From the comparison between Eqs. f!13|) and (JHJ), Ref. [5 
concludes that 

oo 

E 5M n <0 ° ( 15 ) 

n 

(and in fact that it vanishes), to avoid 1/Q 2 terms in the OPE, and, from this, that 

A 2 

5M 2 ~ (-l) n — , n -> oo (16) 
n 

is the slowest possible fall off with n. Equations (fTo^) and @ then lead to Eq. 0. 

Both these conclusions are unfounded. First, Eq. (JT5j) assumes that the sum over n 
in Eq. piijl commutes with the expansion in powers of 1/Q 2 , and hence that the sum 
5^^° converges. 6 In writing Eq. (JTHJ) we assumed that the sum in that equation 
converges, which is indeed the case if SM 2 increases less rapidly than n. But clearly 
this does not imply that &M 2 converges, much less that it vanishes. With the 

5 We ignore as corrections to the parton-model logarithm, and, likewise, the fact that lis and lip 
have non-trivial anomalous dimensions. 

6 We refer to Ref. f° r a detailed discussion on how the sums in Eq. JHJ should be regulated so 
as to be consistent with chiral symmetry. Obviously, the regulator should be taken to infinity before 
taking Q 2 — a physical scale - large [T^] . 
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super convergence relation gone, there is also no ground for the second claim, Eq. (|16|) . 
Asymptotic behavior like that of Eq. may of course be consistent with the absence 
of 1/Q 2 terms in the OPE, but it certainly does not follow from it. 

Second, we also see immediately that demanding that SM 2 be finite leads 
to inconsistent physics, because it would follow that limg2^ 00 U S -p(Q 2 ) is equal to 
AIl(O). The reason this is inconsistent is that chiral symmetry at low energy predicts 
that this constant does not vanish according to Eq. (J7j), while perturbation theory 
predicts that Us-p(Q 2 ) should vanish for Q 2 — > oo (cf. Eq. (|14|)). expressing the fact 
that perturbation theory does not know about chiral symmetry breaking. 

All these problems get resolved by the simple fact that the assumptions made about 
the scalar and pseudo-scalar two-point functions (large- N c , Regge behavior and parton- 
model behavior) can be fully consistent with the OPE, regardless of what happens to 
^2 n SM 2 . We will demonstrate this by explicit example in Sec. 4. 



3 Absence of parity pairs in the 't Hooft model 

Before getting to our example, it is instructive to consider the situation in a simple 
solvable model, two-dimensional QCD at N c = oo [T^. The asymptotic behavior of 
the meson spectrum in this model is known: at large resonance number, the masses 
scale like 

M 2 ~ nA 2 , n -»• oo , (17) 

and resonances of opposite parity alternate: resonances with odd n have the same 
parity, which, in turn, is opposite to that of resonances with even n. 
It follows that the numerator of Eq. ((TJ is given by 

Ml , , - Ml k 

M n+1 -M n = — =±i -p. = , n -> oo , 18 

+ M n+1 + M n 2^ ^ ' 

while the denominator is given by 

M n+2 - M n ~ A= , n -> oo . (19) 



Jn 

In the 't Hooft model, therefore, 

Iim|p(n)| = i (20) 

n— too Z 

While mass differences decrease as a function of n, there is clearly no way in which 
parity pairs can be unambiguously identified. This is reflected in the result for the 
spectral overlap, Eq. ()20j): the value 1/2 is the maximum possible value that p(n) can 
have. This means that asymptotically the pseudo-scalar states lie precisely halfway 
two successive scalar states. 7 In this sense, chiral symmetry is maximally broken by 
all meson resonances in the 't Hooft model, unlike what is claimed in Refs. jHld- 

7 A larger value would imply that parity pairs have been misidentified. 
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We note that the quantity x{ n )-> defined in Eq. (J2J, goes like 

X(n) ~ , n -> oo , (21) 

simply because the denominator of x{ n ) goes like 2-^/nA. The spectrum gets more 
crowded with higher n, and thus naturally the pseudo-scalar states get more squeezed 
between the scalar states. The 't Hooft model thus confirms that x{ n ) is n °t a good 
measure of chiral symmetry restoration. 



4 A counter example in four dimensions 

Returning to four dimensions, we will now construct an example of a mass spectrum 
which satisfies all OPE constraints, but which does not satisfy Eqs. (fT3j) and (|16J) . Our 
model will be that of Sec. 2, with the specific choice 

M 2 Sn>Pn = m% P + nA 2 , n = 1, . . . , (22) 
Fa p = kA 2 MI p , n = 1, . . . , 

where k is a numerical constant, equal to N c /32tt 2 in large-iV c QCD, and we choose 
ms 7^ nip. In the pseudo-scalar sector, there is an additional massless pion pole, with 
residue 2B 2 f 2 (cf. Eq. ©). 
Using 

N / 1 1\ 1 

i im E — — = — ^ , ( 23 ) 

n=l x ' 

in which 7^ is the Euler-Mascheroni constant and i/j(z) is the digamma function 

f°° f e~ l e~ zt \ d 

*<*>=/ * — '^rw, (24) 







we find that 



n 5 _ P (^ = -g 2 ) = (25) 



An(o) -^- + *Q 2 



Q 2 + m 2 \ (Q 2 + m 2 P \ A 2 A 



1> , a -v> 



A 2 7 r V A 2 y Q 2 + m 2 Q 2 + m|, 



An(0) + K(m 2 s - m 2 P ) - (V/ 2 + i«(m| ~m% + A 2 (m 2 s - m 2 p))^j ± 
+k Q(m| - m P ) + U 2 (m* - m%) + l -A\m 2 s - m 2 P )^ ± + O (± 



where we used 



il{z)=log: -l-±^, (26) 

n=l 
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valid for Re z > (B 2n are the Bernoulli numbers). Requiring that this be equal to 
the leading OPE expression, —Q7ias{ipip) 2 /Q 4 , leads to the relations 

16B 2 L S = n(m 2 P - m 2 s ) , (27) 

B 2 f, = ^n(m P -m' s + A 2 (m 2 P -m 2 s )) , 

6na s (^) 2 = lK(2(m%-mi) + 3A 2 (m 4 P -mt)+A\m 2 P -m 2 s )) , 
6 

where we used Eq. (J7J). Clearly, our model wants m 2 P to be larger than m|, which is 
not inconsistent with what is known about the lowest lying resonances in the scalar 
and pseudo-scalar channels. 

This example demonstrates several facts. First, in this model, SM 2 = m| —m 2 P for 
all n = 1, . . . , oo, and therefore Ylm=i $M 2 does not converge. Nevertheless, there is 
no clash with the OPE. Second, the coefficients of the OPE appearing in Eq. ()27|) are 
related to low-energy constants and know about resonances at all scales, and not just 
about those for asymptotically large n. Clearly, if we would include other low-lying 
resonances which are not part of the equally-spaced towers, their masses and residues 
would show up in these sum rules as well. Our example satisfies all the assumptions of 
Ref. but invalidates the claims. In particular, it is straightforward to see that this 
model does not show chiral symmetry restoration: 

/ x • [m 2 P — m% ml + A 2 -171%} . . 

hm p(n) = mm {^^, P ~) (28) 

(assuming that m| < m 2 P < m 2 s + A 2 ). 



5 Discussion 

In this paper, we have restricted ourselves to N c = oo, because in that case all reso- 
nances in a certain channel are stable, and their masses are thus well-defined. Within 
this framework, we have demonstrated that from what we know about the short- 
distance behavior of QCD (in this limit) nothing can be inferred about parity doubling 
in the hadron spectrum. The basic reason for the difficulty comes from the following 
pitfall: even though it is true that perturbation theory tells us that chiral symmetry 
breaking effects die off at high energy, this is true only in the euclidean regime. The 
spectrum lies in the Minkowski regime, and perturbation theory clearly fails badly for 
infinite N c , predicting a smooth continuum where there is in fact an infinite tower 
of Dirac delta functions. The lack of a straightforward connection between the two 
regimes is at the origin of the difficulty. We illustrated the absence of parity doubling 
in detail in the context of a Regge-like model for scalar and pseudo-scalar mesons 
which satisfies both the requirements of large iV c in the Minkowski regime and those of 
leading-order perturbation theory in the euclidean regime. We see no reason why the 
situation would be different in channels with different spin and isospin. 

A further fact illustrated by our example is that the coefficients of the OPE know 
about all resonances in a certain channel, and not only about the asymptotic behavior 
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for large n. This is as one would expect from the OPE, which can be used to connect 
information about hadron spectra as a whole to the behavior of correlation functions 
deep in the euclidean regime. The fact that this connection must work constrains the 
possibilities for the spectrum quite generally. For instance, in Sec. 2 we showed that a 
superconvergence relation of the form of Eq. (J1J cannot hold for large N c , because it 
is in conflict with perturbation theory. 

In the model we considered in Sec. 4, potential candidates for parity doublers can 
still be identified, because the slope in scalar and pseudo-scalar channels is taken to be 
equal, as was done in Ref. 0. This leads to a definite prediction for pin). However, we 
do not really know whether the slope in the scalar and pseudo-scalar channels is equal, 
and one way to generalize the model is to allow for As ^ Ap. Even in this case, a set of 
sum rules generalizing Eq. ()27)1 can still be derived, and thus there is no conflict with 
the OPE. But in this case it clearly makes little sense to try identify potential parity 
partners, let alone define p(n). 

Another important generalization is to go to the world of finite N c . Assuming that 
at finite but large N c the masses of the model of Sec. 4 pick up decay widths T n which 
scale like ~ \/n/N c |Hj, it can be shown |H1 that as a consequence in this model 
the spectral function - Im (n^t) — Up(t)) exhibits a fall off ~ l/£ 2 , for large t. 8 It 
would however simply be incorrect to interpret this fall off in the spectral function as 
a sign of chiral symmetry restoration between the scalar and pseudo-scalar channels. 
It is just the growing of the widths with excitation number which causes this effect, 
by washing out the delta functions in Eq. (J8j). and thus increasing the overlap between 
opposite-parity doublets. The mass spectrum (defined as the location of the real part 
of the masses) is still given by Eq. (|2"2"|) . and no parity alignment takes place. 

While the details of this argument are based on the finite- N c version of our model, 
it is generally expected that broadening of resonances toward larger n takes place 
also in QCD. This smoothens the spectral functions, ameliorating the connection to 
the perturbative regime, and one expects that the difference of the S and P spectral 
functions will fall off bounded by some power of t also in the real world. However, 
our model demonstrates that this does not imply that chiral symmetry restoration has 
to take place. The slowest possible pattern of chiral symmetry restoration compatible 
with the OPE is no chiral symmetry restoration at all. 
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